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Abstract. We investigate the interplay of quark and meson degrees of freedom in a physical state repre- 
senting a near-threshold resonance for the case of multiple continuum channels. The aim is to demonstrate 
the full complexity of possible near-threshold phenomena. It turns out that those are especially rich, if both 
quark and meson dynamics generate simultaneously weakly coupled near-threshold poles in the S-matrix. 
We study the properties of this scenario in detail, such as t-matrix and production amplitude zeros, as well 
as various effects of the continuum channels interplay. 
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1 Introduction and Main result 

A convenient parametrisation of an 5- wave near-threshold 
amplitude was suggested in Ref. [Tj . Consider an unstable 
particle coupled to the hadronic channel, open at E = 0, 
with the coupling constant gj. The t- matrix in the effec- 
tive range (Flatte) approximation can be written in the 
form: 



t(k,k,E) 



1 



9f 



1 



fJJ 



8n 2 fi E - E f + \g s k 8n 2 fiV(E) : 



(1) 



where 



k(E) = y / 2jlE0(E) + iy/-2nEG{-E), 



and fi is the reduced mass in the hadronic channel. Nor- 
malisation of the t-matrix is such that the scattering am- 
plitude is given by 

F(E) = -4n 2 fj,t(k,k,E). 

In fact, the Flatte form of the scattering amplitude is noth- 
ing but the Breit-Wigner amplitude with the momentum 
dependence of the elastic width taken into account ex- 
plicitly. However, as argued in Ref. [2], direct interactions 
between mesonic channels could lead to quite a peculiar 
form of the line shapes, drastically different from the ones 
driven by the simple Flatte form (fTJ). 

The simple form (fTJ) is valid in the near-threshold re- 
gion, that is, for 

\E\ « \A\, (2) 
where A denotes the distance to the next threshold, 



with M t hi being the i-th threshold. However, in practise, 
analysis of experimental data requires knowledge of the 
amplitude in a wider range, covering a few thresholds. In 
this case the simple form (JTJ has to be modified. Con- 
sider two isospin-related channels. A naive generalisation 
of Eq. (JTJ to this case would reacQ: 



V(E) = E-E f + -g f (k 1 + k 2 ), 



(3) 



where k% and k 2 stand for the momenta related to the 
corresponding thresholds: 



ki(E) = y/2jiEO(E) + iy/-2nEO(-E), 
k 2 (E) = y/2 f iE - A)0(E - A) 

+ i^2 l x(A-E)9(A-E), 



(4) 



and further generalisation to the multi-channel case is triv- 
ial. However, if direct interactions between mesonic chan- 
nels are present in the system, the simple form (J3J) may 
become inappropriate and thus it may be subject to gener- 
alisation, quite analogous to the one-channel case already 
discussed in Ref. [2j. Usually the analysis of experimen- 
tal data based on the multi-channel Flatte parametrisa- 
tion of the near-threshold observables still employs simple 
expressions similar to Eq. ([3]). with all effects of the di- 
rect interaction between mesonic channels neglected or, 
more generally, parametrised into the effective coupling 
constants — see, for example, Refs. [3J|3J[S]. It is argued 



Notice that, 



assuming isospin conservation, we set = 

gy — gf and that here and in what follows the constant gj 
differs by a factor of two from the one introduced in Eq. (fTJ) . 



A = Mi 



U,0 



M, 



thi ■ 



,(2) 
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in Ref. [5], however, that the explicit inclusion of such ef- 
fects could lead to line shapes which are even more unusual 
than those obtained in the one-channel case discussed in 
Ref. [2] . The aim of the present paper is to study system- 
atically the role of the direct interaction between mesonic 
channels. In particular, we investigate under which con- 
ditions the effect of direct interactions on the line shape 
cannot be parametrised into the couplings and has to be 
taken into account explicitly. Throughout this paper we 
assume the constituents of continuum channels as stable 
particles. Although a finite width of a constituent can have 
a significant impact on line shapes, as already discussed, 
for example, in Refs. [IKE, this effect can be treated 
independently from the interplay of compact states and 
continuum channels, which is the focus here. 

As the main quilitative result of our research we find 
that in most cases one should expect the multi-channel 
Flatte parametrisation of Eq. ^ to describe properly the 
line shape of near-threshold resonances, however, under 
certain conditions (especially for fine tuning of parame- 
ters) more complicated line shapes can emerge. We thus 
propose the following procedure for the data analysis: (i) 
if the form of Eq. ([3]) provides enough structure to de- 
scribe new-threshold data, more complicated parametri- 
sations cannot be used in a mo del- independent way — the 
near-threshold data are simply not sensitive enough to dis- 
criminate between different parameter sets which provide 
identical near-threshold behaviour and produce different 
structures only far beyond the near-threshold region; (ii) 
if, on the contrary, line shapes in the near-threshold region 
demonstrate irregularities, more complicated structures, 
as given below, should be used. In this case more than 
two poles are located near the relevant threshold and a 
non-trivial interplay of various effects is observed. If this 
would ever be relevant in practice will depend strongly 
on how well the background is understood for the partic- 
ular measurement; (iii) finally, the mulit-channel formal- 
ism developed below may be used to study the resolu- 
tions/statistics necessary for the data to actually detect 
the possible interplay of compact states and continuum 
channels. 

The paper is organised as follows. First, in Sec. [5] we 
derive general expressions for the t-matrix and the wave 
function (w.f.) of a near-threshold resonance which appear 
as a result of the interplay of a bare quark state and mul- 
tiple continuum states. This section is quite technical. It 
forms the foundation of our work but can be skipped by 
those who are more interested in examples of phenomeno- 
logical applications of the formalism developed. Those are 
given for the two-channel, near-threshold case in Sec. |31 
which is then studied in detail from the point of view of 
the effect of the direct interaction between mesonic chan- 
nels. The formulae derived are exemplified further by the 
case of the A (3872) charmonium in Sec. |U In particular, 
we argue that description of the experimental data for the 
A (3872) charmonium in the near-threshold region does 
not require severe modifications as compared to the naive 
approach and can be safely performed with the help of the 
simple formula Finally, in the concluding part of the 



paper, we discuss various practical aspects of the use of 
the formulae derived in the data analysis. 



2 The multi-channel problem 

2.1 Essentials of the formalism 

We study a physical state (for the sake of convenience 
referred to as A) which is a mixture of a bare qq state 
(labelled below as 0) and multiple two-body components 
(labelled as i = 1, 2, . . .) and represent its w.f. as: 

/ c|Vb) 



I*} 



V 



Xi(p)\MnM 12 ) 
X2(P)|M 21 M 22 ) 



(5) 



with p being the relative momentum in the corresponding 
meson pair {MiiM&} • 

We allow for direct interactions between mesonic chan- 
nels i and j encoded in the potentials Vij(p,p') (including 
diagonal terms with i = j). Transition between the i-th 
mesonic channel and the quark state is caused by the tran- 
sition potential fi(p). Then w.f. ([5]) obeys a Schrodinger- 
like equation: 

H\X) = M\X), (6) 
with the Hamiltonian 



H = 

Here 
and 

H hi \p,P'] = 



(7) 



f H A[p] f 2 [p] —\ 
h[p] H hl [p,p'] V l2 [p,p'] •• 

Up] v 2 i[p,p'} H h2 [p, P '}-- 

v ) 

H \ip ) = Mo|Vo) 



m a + m l2 + ^— ) <5 (3) (p-p') + % [p, p'] 



where mu and m^ are the meson masses in the i-th chan- 
nel, with the reduced mass 

m ix m i2 

(M = ; ■ 

mil + m i2 

We use hats for the operators and square brackets for their 
arguments in order to distinguish between operators and 
the corresponding c-number functions. In particular: 

6[p}c=0(p)c, 

0[P] X (P) = J 0(p) X (p)d 3 p, 

d[p,p']x(p) = J o(p,p') x (p')d 3 p'- 

In what follows we define the energy E relative to the 
lowest threshold: 

M = mn + mi2 + E, M = mu + m 12 + E - 

Then the mass parameters 

A = (ma +m i2 ) - (mn +mi 2 ), (8) 

define splittings between mesonic thresholds. 
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2.2 t-matrix should be stressed that the procedure used to arrive at the 

solution (ITT1) - (fT4|) automatically ensures unitarity, since it 
The system of coupled Lippmann-Schwinger equations for is based on solving a Schrodinger-type equation, 
the various components of the off-shell i-matrix reads: Finally, for further convenience, we define quantities 

Ri and R[ such that 

= / fk(q)S k (q)t k0 (q,E)d 3 q, , 

k J / fi(q)S I (q)d 3 q = /* (J2, + Hi), 



t i0 (p,E) = fi{p) 



V l k(p,q)S k (q)t k0 (q,E)d 3 



t 0i (p,E) = f i (p)-J2 I f k (q)S k (q)t kl (q lP ,E)d 3 q, 



where f oi = /;(0) and 

Ii = 47T 2 /X i fc i , 



fi(p)t 0j (pf,E) 



tij(p,p', E) = Vij(p,p') l — 3 — — — — with the momenta h defined similarly to Eq. Q. that is 

Hjq — hi — ?u 

-J2fVik(p,q)S k (q)tk j (q,p',E)d a q, ^(E) = y/2iu(E - Ai)0{E - A t ) (19) 

k " (9) +iy/2n i (A i -E)e(A i -E). 

where ^ 

^\P) — p2 /(2fii) — E + A~i — iO 2.3 Bound states and continuum spectrum 

is the propagators of the meson pair in the i-th channel. The Schrodinger-like equation © admits both bound- 
Solution of system © can be found in the form (see also state solutions and continuum spectrum. If a bound state 
Ref. my- is present in the spectrum of the Hamiltonian (JT]), the 

corresponding solution of Eq. © takes the form: 

t <F)= (E-Eo)g(E) 

° 0{ ' E-Eo+GiE) 1 1 ; c B = cos6, X iB =Si(p)&(p) cos 8, (20) 

toi(p,E) = — — (j>i (p), (12) where the hadronic Green's function Si(p) and the func- 

E - E + 9(E) tion are defined in Eqs ^ (jjoj and (Qi^ respectively. 

t o(p E) = E — Eq The value of the angle follows from the normalisation 

1 ' E — E Q + Q(E) 1 ' condition for the w.f. ^ which, for the bound state solu- 



M^> = «<^ + ^^> M4, 



tion, reads: 



where 



E / d 3 p\ XlB {p)\ 2 = 1, 



MP) = fi(P)-J2[ t yk(P^)S k (q)f k (q)d 3 q, (15) S ° ^ /" 

k J tan 2 # = E / S i(Pm(P)<rP- 

4>i{p') = fi(p) - E/ S k(q)fk{q)t kl (q,p')d 3 q, (16) The Weinberg ^factor (see Ref. 0), which measures 

fc the probability to find the compact component in the 

g^g^ _ \ / f 2 ( y q')S (q)d 3 q (17) physical w.f. and therefore encodes information on the na- 

. J % 1 ture of the state, is then: 

- E/ Mk)S t (k)t^(k, q )S 3 (q)f 3 (q)d 3 kd 3 q, Z = \c B \ 2 = cos 2 6. (21) 

The i-matrix of Eqs. ([TTjl - fTl)) develops a pole at the 

and the components of the t-matrix t v (p, p') by definition bound-state energy e (that is, at Af B = m u + m 12 - e), if 
satisfy the last equation in system flS} with hip) = 0, that p rv \ fo*^ 

lb, 

„ where the loop function Q (E) is defined in Eq. (fTT) . Alter- 

tYj(P>P') = Vij(P,P') ~E / ^fcGP>g)£*(«)*&(9>P0A. natively, Eq. {25) can be obtained directly from Eqs. © 

k J and ((2U)) . The generalisation to the case of multiple bound 

(18) states is trivial. 

and thus the t- matrix t v (p,p') absorbs all the details The Hamiltonian of Eq. ([7]) possesses also a continuum 

of the direct interactions between mesonic channels. It spectrum. The solution of the Schrodinger-like equation 
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^ with free asymptotics in the hadronic channel i takes 
the form: 

X$h(p) = <M (3) (P - ki) - S^t^kuE), (23) 



t 0i (ki,E) 
En-E ' 



(24) 



where ki — y/2fii{E — Ai). Equation (|2~4"|) allows one to 
calculate the spectral density w(E), 



;(£)=£> ifei 0(£-A) / |c£(£) 



dot 



(25) 



which gives the probability to find the bare state in the 
continuum w.f. with the given energy E [5] . Then the over- 
all normalisation condition reads: 



/>oo 

Z+ w{E)dE = 1, 
Jo 



where Z is defined in Eq. (f2~T|) . In case of existence of 
several bound states, the latter equation is obviously gen- 
eralised as 



/•OO 

VZ„+ / w{E)dE=l, 
Jo 



where the corresponding iJ-factors are given by expres- 
sions similar to Eq. (|2~T1) . 



3 Two-chanel case 

In this chapter we present the simplest nontrivial appli- 
cation of the general formulae (fTTjl -tfH I derived above, 
namely, we consider the two-channel solution and perform 
its low-energy reduction. We have A\ = and A 2 = A > 
0. In what follows, assuming smallness of A as compared 
to the masses, we set fii = fi 2 = M- 



3.1 i-matrix, bound states, and Z-factors 

In the scattering length approximation the 2x2 direct in- 
teraction i-matrix t v (p,p') (see Eq. (JTSJ) for its general 
definition) can be parametrised in the form: 



t v = 1 (\{R s + R t ) + ih \(Rt-R s 



Det 



\{R t -R s ) \{R s + R t )+ih 



(26) 



where 



Det = (R s R t - hh) + ^{Rs + Rt)(h + h), (27) 

and the quantities R s and Rt can be related to the (in- 
verse) scattering lengths in the singlet and triplet chan- 
nels, respectively: 



R s = 4ir 2 fij s , R t = 47r 2 /x7 t . 



(28) 



Note, those are not the physical scattering lengths, but 
only the scattering lengths that would have been realised, 



if only the potential scattering had been present; in general 
the physical scattering lengths are combinations of 7 S , 74 
and the pole parameters. In what follows we stick to 7 S 
and 7t only as to convenient quantities parametrising the 
strength of the direct interaction i-matrix t . 

In what follows, we assume that the quark state is an 
isosinglet, setting 

R,\ = R'2 = R, R[ = R' 2 = R', foi = fo2 = fo/ V%, 

and introduce three physical parameters, Ef (Flatte en- 
ergy), Eq (zero of the singlet i-matrix — see below as well 
as Ref. 0), and gf (Flatte coupling constant) instead of 
the three bare parameters Eq, R, and R': 



f 2 a 
Ef — E — —(RR S — R ), 

Ec = Eq — /q(R s + R — 2i?'), 



9f 



8tt> 
R 2 



fo(R s - R' 



(29) 
(30) 
(31) 



Then, after straightforward calculations, one can find 
for the various components of the t-matrix: 



t 



11 



Ui = t 



1 -f s (E-E f ) + (E-E c )ht + 2ik 2 ) 
8ir 2 n D(E) 

1 lt {E - E c ) - ls{E - E f ) 



8n 2 fi D{E) 
1 ls {E-Ef) + {E-E c ){-it + 2ifci) 



8tT 2 /1 



D(E) 



(32) 
(33) 
(34) 



where 



D(E) = 7s ( 7i + 2 (fci+fc) ) [E-Ef) 



k 1 k 2 -- lt {k l + k 2 )){E-E c ). (35) 



Expressions (|32|) — (f5^|) coincide with those derived in 
Ref. [B] if the parameters 70 , 71 , and g used in Ref. [BJ are 
mapped to our parameters as 



1 



7o= 7s, 7i=7t) 



Notice also that a different definition of the isospin states 
is used in Ref. [6] (see Eq. (|49| below), so that the t- 
matrix components in and £22 appear in Ref. [BJ with the 
opposite signs as compared to our Eqs. (|3"2")l and (|34p. 

It is instructive to construct the isospin singlet and 
triplet combinations, which read: 



t„ = 



tt = 



tst 



1 (E - E c ) (2 7t + ijjei + fc 2 )) 
8ir 2 n D(E) 

1 2 ls (E-Ef)+i(k 1 + k 2 )(E-E c ) 
8ttV D(E) 

1 i(k 2 - h)(E - Ec) 
8ttV D(E) 



(36) 
(37) 
(38) 
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where, similarly to the single- channel case [2], the i-matrix 
t s possesses a zero at E = Ec- In addition, expressions 
(j3"rj)) - ([3"51) contain explicitly the leading isospin violation: 
k\ 7^ k 2 allows for transitions between isospin singlet and 
triplet states. 

For the bound-state solution of the Schrodinger-likc 
equation ([6]) the Z-factor can be found explicitly in the 
form: 

Z=- lr-^ = - (39) 



l + 0'(-e) l + tan 2 6»' 



with 



(40) 



tan'g=^ c + e) 
k,\k 2 

Kl(jt - Kx) 2 + K 2 (7t ~ K 2) 2 

(ki + k 2 - 2j t ) (j s jt + kik 2 - \{ls + 7t)( K i + K2)) ' 



where 
case. 



Im(fcj(— e)) — see Eq. (|19p for the general 



3.2 Limiting cases 

In this chapter we check various limiting cases of the low- 
energy two-channel formulae derived above in this Section. 

We start from the single-channel case. In order to ap- 
proach the single-channel limit one is to set A = every- 
where which, in particular, implies the substitution 

ki(E) — k 2 (E) = k(E). 

Then one readily finds for the singlet t-matrix and for 
the Z-factor: 



1 - Z 

z 



tan 2 6 



E - Ec 



E c + s 



(41) 



2s 



■Is 



respectively, which coincide with the corresponding ex- 
pressions derived in Ref. [2] , after an obvious identification 

Is = IV- 

The zero Ec affects the dynamics of the system, if 
it appears in the near-threshold region. For the coupled 
channel situation at hand this happens for \Ec\ ^ A 
c.f. the discussion in Ref. [2] for the single-channel case. 
The value of Ec is, by construction, controlled by the 
singlet inverse scattering length 7 S , since it follows from 
Eqs. P9 l) -(|5T ]) that 



Er 



E 



1 



/ 



;9fls 



(42) 



Thus the zero Eq is generated far away from the near- 
threshold region of interest, \Ec\ ~> A, if 



. A 

hs » — . 

91 



(43) 



For this estimate we assumed \Ej\ < A. We therefore con- 
clude that, if Eq. (|4"5)) holds, then 7 S might be assumed to 



be infinite or, correspondingly, that the singlet scattering 
length is negligibly small. It should be stressed however 
that this does not imply the absence of meson-meson in- 
teraction at all in the singlet channel. Indeed, as explained 
before (see the discussion below Eq. ([2"5])). the scattering 
lengths 7 S and 7* are only residual scattering lengths. Sim- 
ilarly to the single-channel case we expect the line shapes 
and the interpretation thereof to simplify tremendously 
for large values of the Ec, however, as we shall see, ad- 
ditional structures can be introduced through the triplet 
channel. 

Indeed, the momenta k\ and k 2 enter expressions for 
the i-matrix not in the form of a naive sum ki + k 2 , 
but in an entangled way, including the product k\k 2 (see 
also Ref. [6] for the discussion of the related effect on 
the line shapes for the X(3872)). This entanglement of 
mesonic channels is governed by the triplet inverse scat- 
tering lengths 7t and it becomes strong for |7 t | < \J [iA. 
Therefore, the largest limit is achieved for 

\lt\ » VT^ (44) 

in the general two-channel formulae. 

We therefore identify four possible cases: 

— Case (i): |*y s j — > 00 and |7t| — > 00. 

— Case (ii): small 7^ and |7t| — > 00. 

— Case (hi): |7 S | — > 00 and small 74. 

— Case (iv): both 7 S and 74 are small. 

To approach Case (i) above one is to assume both con- 
ditions (|4U| and ([4"4"]) to hold simultaneously. In this case 
one arrives readily at the two-channel Flatte denominator 
which contains trivially the sum of the two channels, that 
is, at Eq. ([3]). For the Weinberg Z factor we find: 



Z= 1 



fj.g f (K 1 +K 2 ) 
4ki« 2 



with the standard pattern (c.f. Ref. [TU]), namely for 3/ — > 
we have Z 1, which corresponds to a compact state, 
while for gf — > 00 we have Z — >• and the state is to be 
interpreted as dominantly moleculafl 

Case (ii) implies that, while entanglement of the meso- 
nic channels still remains irrelevant, the zero Ec starts to 
approach the near-threshold region. Then, in this limit, 
one arrives immediately at the trivial generalisation of the 
single-channel expressions found in Ref. 2 . In particular, 
the Flatte denominator takes the form: 



V F {E) = E — Ef ~ 



(E - Eff 



E-E, 



c 



+ -gfih + k 2 ), (45) 



which simply amounts to the substitution k — > \{k\ + 
k 2 ) in the corresponding single-channel formula. Then the 
entire analysis performed in Ref. [2] with respect to the 



It is important to stress that from the Flatte parametri- 
sation alone it is not possible to reliably extract gf, if it is 
large, due to a scale invariance of the parametrisation in that 
limit pi]- 
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role played by the zero Eq for the dynamics of the system 
applies directly to the two-channel case at hand. In order 
to keep in Eq. (|45[) the first non-vanishing correction in 
the l/jt expansion one is to substitute: 



fci 



fci + fc 2 



27t 



In Case (iii) , the zero Eq plays no role in the dynamics 
of the system, the new effect as compared to the simplest 
Case (i) being the entanglement of the mesonic channels 
governed by the triplet (inverse) scattering length 7t . The 
corresponding expressions for various components of the 
i-matrix read then: 



t, = 



9f 1 
8% 2 fiV (£)' 

E-E f + {g f {k 1 + k 2 ) 1 



U 47r2/i( 7t + £(fci + fc a )) V{EY 
, igf k 2 - ki 1_ 

st ~ i6^ 2 /i 7t + i(fc 1 + fc 2 )i?(i;)' 



(46) 
(47) 

(48) 



where 



The effect of the channels entanglement for the A (3872) 
resonance was discussed in Ref . [5] . As a result new struc- 
tures in the corresponding line shapes were observed es- 
pecially far away from the threshold, at E A. In this 
paper we choose to stay within the range of applicability 
of the effective-range expansion used to derive all formulae 
above, that is in the near-threshold region \E\ < A. 

Finally, in Case (iv) one cannot simplify the general 
formulae derived in the previous chapter and is forced 
to face the full complexity of the interplay of quark and 
mesonic degrees of freedom as well as that of the mesonic 
channels entanglement. 

One should have in mind that, in reality, inelastic chan- 
nels with far-away thresholds also exist, and could be dealt 
with the help of the general multi-channel formalism of 
Sec. [5] However, these remote thresholds do not produce 
extra energy (momentum)-dependent contributions to the 
near-threshold amplitude (though their presence in the 
system may modify the line shapes appreciably — see, for 
example, Ref. [4]). The net effect of such thresholds can be 
taken into account by including imaginary parts of inverse 
scattering lengths 7 S and 7 t (if the inelastic channels are 
coupled to hadronic components of the w.f.) and Flatte 
energy Ef (if the inelastic channels are coupled to the 
quark component). Since the present research is aimed 
solely at the discussion of multi-channel dynamics in a 
near-threshold resonance, we regard the investigations of 
remote thresholds as lying beyond the scope of this paper. 



4 An illustrative example 



which incorporated effects due to quark and multi-channel 
mesonic dynamics. For the two-channel case we performed 
the low-energy reduction of the solution found and arrived 
at rather simple expressions. We then pinpointed parame- 
ters which governed the dynamics of the system, in partic- 
ular, the appearance of the zero Ec in the near-threshold 
region and the entanglement of mesonic channels. This 
finalised a purely theoretical analysis of the interplay of 
quark and meson degrees of freedom in a near-threshold 
resonance in the multi-channel case. In this chapter, we 
proceed with phenomenological applications of the gen- 
eral formulae derived above. In particular, we study line 
shapes of the A production under various assumptions 
concerning the production mechanism. Keeping in mind 
that the charmonium state A (3872) provides a paradig- 
matic example of a two-channel situation, we use it as an 
anchor which helps us to fix several sets of parameters. 
We caution the readers however of a naive identification 
of these sets of parameters, used solely for illustrative pur- 
poses, with real fits to the experimental data for the state 
A (3872). As was explained before, the analysis of the ex- 
perimental data for the A (3872) lies beyond the scope of 
the present paper, since the aim of the latter is develop- 
ing and improving a phenomenological tool which could 
be used then for the analysis of experimenal data on vari- 
ous near-threshold resonances. Three more comments are 
in order here before we proceed. 

First, we fix the quantum numbers of the A to be 1 ++ 
which is compatible with the A (3872) state. Let us remind 
the readers that the problem of the quantum numbers for 
the A(3872) is not fully resolved yet: while the analyses 
of the J/ipTr + ir~ decay mode of the A (3872) yields ei- 
ther 1 ++ or 2 h quantum numbers [T2], the recent anal- 
ysis of the J J 'ifm + tt~ 7T° mode seems to favour the 2 h 
assignment [13], though the 1 ++ option is not excluded. 
As shown in Ref. [2], the A (3872) cannot be a naive 
cc 2 ^ state and, were the 2 h quantum numbers con- 
firmed, very exotic explanations for the A (3872) would 
have to be invoked. In the absence of such a confirma- 
tion we stick to the most conventional 1 ++ assignment for 
the A (3872). With the given quantum numbers, the low- 
est open-charm channels are Z?L>i3i n the relative S'-wave, 
and the X (3872) is known to be located very close to the 
D°D*° threshold. Identification of the channels is obvi- 
ous in this case: channel 1 is associated with the neutral 
component D°D*° while channel 2 — with the charged 
component D + D*~. 

As the second comment we would like to note that, 
for the definition of the C-even(odd) state, we stick to the 
convention of Ref. [TS] for the C-parity operator C and 
C-parity eigenstates: 

P(V) = CP(V), C\PV±PV) = ±\PV±PV), (49) 

with P and V denoting the pseudoscalar and the vector, 
respectively. Further details of the ambiguity in the defi- 
nition of C-eigenvalue states can be found in Ref. [T5] . 



In previous chapters we derived general expressions for 3 An obvious shorthand notation is used here and in what 
the w.f.'s and i-matrices for a near-threshold resonance follows: DD* = -^(DD* + DD*) (see Eq. ((49)0 . 
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The third comment concerns the stability of the X con- 
stituents Mix, M 12 , M 2 i, and M 22 (see Eq. ([5]) above). 
In order to get an insight into this problem we resort 
again to the example of the charmonium A (3872). The 
D* mesons are unstable, and the question is whether the 
account for the D* finite width is important. On one hand, 
if the A (3872) is a bound state and thus its pole is located 
on the first sheet below the D°D*° threshold, then finite 
width effects produce a spectacular below-threshold peak 
in the D°D°ir° mass distribution — see Refs. [4,5 . As 
argued in Ref. [1], observation or non-observation of this 
below-threshold peak is a potentially powerful tool to dis- 
criminate between various models for the A(3872). On the 
other hand, the existing data analyses does not allow one 
to exploit this tool. As stressed in Refs. [5], the analyses 
procedure is biased by the assumption that all D°D°ir° 
events come from the D°D*° mode, so that events from 
below the D°D*° threshold are move above it. In addi- 
tion, the existing experimental resolution is too coarse to 
investigate in detail the narrow region around the D°D* a 
threshold. Finally, the finite width effects are known to 
be negligible above threshold for the energies larger than 
the width of the D*°, that is, larger than approximately 
100 keV (see the discussions in Refs. [H|5j[Z])- Therefore in 
what follows we neglect these effects. 

We stick to the neutral mesons in the final state and 
consider production of the X through the hadronic as well 
as the quark components. The corresponding amplitudes 
read: 



M-hi — J'h! — Fhx 

1 



-To 



E-E 



S 1 (p)t n (p,k,E)d 3 P , (50) 
S 2 (p)t 21 (p,k,E)d 3 p, (51) 
toi(E), (52) 



where E = k 2 /(2fi), and T's denote point-like sources 
which incorporate the details of the short-ranged dynam- 
ics of the X production. In what follows we neglect Born 
terms in the production through the hadronic components 
(the first term on the r.h.s. of Eqs. ((511)) and (|5i~j0 . This 
is justified if the t-matrix elements have a near-threshold 
pole (a detailed discussion can be found in Ref. 0). 

Then, if one neglects the interference between three 
possible production mechanisms, it is straightforward to 
find the corresponding branchings: 



dBr « = 0(E) Bo 



dE 
dBr hl 
dE 



dBr h . 2 
dE 



6(E) 



6(E) 



D(E)\i 
Biy/E 



7 S 2 \j t + ik 2 \ 2 . 



D(E)\i 

+ (it 

b 2 Ve 

D(EW 



Is (E - E f ) 
+ 2ik 2 ) (e - 
Is (E - E f ) 



(53) 
(54) 



E f + 2^/7* 



(55) 



-It 



E-E f + -g f j s 





7 8) MeV 


7t, MeV 


E f , MeV 


Line 


Z 


(i) 


±00 


±oo 


-10.47 


solid 


0.30 


(ii) 


-30 


±oo 


-3.22 


dashed 


0.85 


(iii) 


±00 


-30 


-7.77 


dashed 
-dotted 


0.19 


(iv) 


-30 


-30 


-2.97 


dotted 


0.67 



Table 1. Scattering length parameters, the Flatte parameter 
Ef, and the Z-factor for the bound- state at e = —0.5 MeV for 
Cases (i)-(iv). 



where the denominator D(E) is given by the following 
expression (see Eqs. (|3"S")) and (H3)): 



D(E)= ls ^ /t + -(k 1 +k 2 )j (E-Ef) (56) 

hk 2 - l - lt (ki + k 2 ) j (e - Ef + X -g nt 

The coefficients Bq, B\, and B 2 are proportional to the 
short-range factors |.Fg| 2 , |^7ii| 2 , and I-FmI 2 and, in the 
case of production via hadronic components, include the 
effects of renormalisation of loop functions in Eqs. (f5T)f 
and (|5ip . The formulae for the production rates can be 
obtained from the general ones given in Ref. j6], if one 
assumes the dominance of certain production mechanism, 
and neglects the finite width effects. 

In order to exemplify the line shapes given by Eqs. 
(f5"3"]) - (|5"5"l) we fix parameters as 



fx = 966.5 MeV, A = 8.1 MeV, g f = 0.25 (57) 

and consider Cases (i)-(iv) introduced before. Other pa- 
rameters are quoted in Table [T] Notice that, for infinite 
inverse scattering lengths 7 S and jt their signs obviously 
play no role. Finite 7 S and 7* are chosen to be small and 
negative, providing in such a way a strong attraction in 
the system which is, however, not strong enough to gen- 
erate a bound state by itself. Extra attraction is provided 
by the coupling to the quark state and, for each case, the 
Flatte parameter Ef is tuned to provide a bound state at 
E B = -e= -0.5 MeV. 

In Figs. [TJ2] we plot the line shapes for the X pro- 
duction through the quark component (see Eq. (|53[) and 
Fig. [I| and both hadronic components (see Eq. (I5"4"|) and 
Fig. [2] as well as Eq. (|55|) and Fig. [3l respectively) for 
all four above Cases (i)-(iv). For each individual curve in 
Figs. 0312] the integral over the near-threshold region (cho- 
sen to be from to 10 MeV) equals to unity (in the units 
of MeV" 1 ) which fixes the overall factors Bo, B\, and B 2 . 



Finally, in Tableland Fig. 21 we show pole positions 
for all four Cases (i)-(iv). To this end we introduce the 
plane of the complex variable ui — see Appendix [X] for 
the definition of the variable lo and for the mapping on 
the Riemann sheets. 
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± 0.63-U.06 
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(-0.96) 




(1.94=f il.54) 
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il.28 


±0.70+i0.19 


-i0.99 




0.3 


. \ \ 

\ 




(-0.5) 


(8.26^11-42) 


(-6.6xl0~ 4 ) 








(iv) 


il.28 


±0.25+i0.73 


±0.9540.31 


±0.50-il.08 


0.2 






(-0.5) 


(0.42=Fil-32) 


(7.73±i0.02) 
-i0.93 


(1.24=pil.l0) 


0.1 










(-0.05) 





E [MeV] 



Fig. 1. Production rate through the quark component (|53[) . 
Cases (i)-(iv) are given by the solid, dashed, dashed-dotted, 
and dotted lines, respectively. 



dBr^/dE [MeV ] 
0.25 




E [MeV] 



Fig. 2. Production rate through the first hadronic component 
(I54|) . Cases (i)-(iv) are given by the solid, dashed, dashed- 
dotted, and dotted lines, respectively. 



dBr,,,/dE [MeV -1 ] 
0.4 




E [MeV] 



Fig. 3. Production rate through the second hadronic compo- 
nent (|55p . Cases (i)-(iv) are given by the solid, dashed, dashed- 
dotted, and dotted lines, respectively. 



A comment on the factorisation approximation used 
before in the production rates (I53|) -(|55 p is in order here. 
For the X(3872), this factorisation approximation for the 
D°D*° production via the hadronic component implies 
that in the charged B-meson decays, the X is produced 
through the neutral component 1 while, in the neutral 



Table 2. Positions of near-threshold poles in cj-plane (in E- 
plane in MeV) for Cases (i)-(iv). The column labels I-IV mark 
the corresponding Riemann sheets. See also Fig. U for the 
graphical presentation of the poles. Notice that some poles 
present in Fig. [3] are omitted in the table since the correspond- 
ing energies lie far beyond the near-threshold region of interest. 



S-meson decays, the X is produced through the charged 
component 2. The X production through the quark com- 
ponent is, obviously, the same for charged and neutral 
JS-meson decays. 

One could have naively concluded that, with data on 
the X(3872) production in both charged and neutral B- 
meson decays available, it is possible to investigate the 
full complexity of quark and hadronic dynamics respon- 
sible for the nature of the X resonance. In reality, how- 
ever, both quark and hadronic mechanisms contribute and 
interfere and, with a suitable choice of the production 
parameters .F's, the resulting near-threshold line shapes 
could be smooth enough even in the presence of compet- 
ing quark and meson dynamics — see the examples given 
in Ref. [5]. We emphasise once again that the line shapes 
depicted in Figs. [T][3l where each figure originates from a 
particular production mechanisms and thus interferences 
are neglected, should not be viewed as realistic and suit- 
able for the analysis of experimental data. With the sets of 
parameters and assumptions used we more aim to indicate 
the full complexity of near-threshold phenomena coming 
from the interplay of the quark and meson degrees of free- 
dom as well as from the direct interaction between mesonic 
channels in an X-like resonance. Notice that in the par- 
ticular case of the X(3872), the production line shapes 
do not exhibit irregularities in the near-threshold region 
(up to possible existence of a bound-state peak which, 
as discussed above, is not observable due to the current 
situation with data analyses and the experimental resolu- 
tion). The (3872) line shapes can be described therefore 
by simple Flatte formulae, which reduce to those of the 
scattering length approximation, when the effective cou- 
pling constant is large (see Refs. [IT] for the details). Such 
a scenario is exemplified by our Case (i) above: there is 
only one near-threshold pole corresponding to the bound 
state, which is predominantly molecular. On the contrary, 
our Cases (ii)-(iv) yield several near-threshold poles in 
addition to the bound-state one, and the line shapes are 
strongly distorted either by the i-matrix zero or by the 
channel entanglement — see Tableland Fig. [U Such ex- 
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Fig. 4. Near-threshold poles in the cj-plane for Cases (i)-(iv). 



-0.5 



Case (iv) 



otic scenarios seem to be largely excluded by the existing 
data on the X(3872). 



5 Summary 

In the mult i- channel situation the dynamics of the reso- 
nance is affected by both the interplay of the quark and 
meson dynamics as well as by the direct interactions be- 
tween mesonic channels. As seen from the examples given 
above both these effects may give rise to quite peculiar 
properties of the resonance line shapes. Thus one con- 
cludes that, if the experimental data exhibit such prop- 
erties, the resonance is generated by this complicated in- 
terplay. The converse is not necessarily true: the interfer- 
ence between different production mechanisms is possible, 
which could tame the resulting line shapes. This, however, 
requires some fine tuning for the interference between dif- 
ferent production mechanisms. 

Our studies allow one to outline a sensible procedure 
to investigate the properties of near-threshold resonances: 

1. If data exhibit irregular behaviour in the near-thresh- 
old region, the latter may be attributed to the inter- 
play of the quark and meson degrees of freedom, and 
formulae for the general Case (iv) above or, possibly, 
for one of its limits, as given by Case (ii) or Case (iii), 
should be used in the analysis. 

2. If data in the near-threshold region does not exhibit 
any irregularities, simple Flatte approximation — see 
Eq. ((3]) — is to be used. Indeed, although fine tuning 
of parameters in the general near-threshold formulae 
is possible in this case, which tames the resulting line 
shape, staying close to the threshold, it is not possi- 
ble, as a matter of principle, to discriminate between 
different parameter sets which produce identical near- 
threshold behaviour. Then, in order to get insight into 
the nature of the resonance one is to analyse data ei- 
ther far from the threshold region or data for other 
production channels for the given resonance. In the 
former case, staying far from the threshold, one has to 
face the full complexity of the multi-channel formal- 
ism described in Sec. [5] above. No model-independent 
analysis is possible in this case and one has to rely 



upon model-dependent assumptions in order to arrive 
at simpler expressions. 
3. If no data are available far from the threshold, the 
general formalism of Sec. [2] can be used in order to 
study to what extend (statistics, resolution, binning 
procedure, and so on) the data would need to improve 
in order to get sensitive to the structures potentially 
present. 

No complicated structures are observed in the near- 
threshold region in the existing data on the X(3872) pro- 
duction in S-meson decays. So, employing the analysis 
procedure outlined above, one concludes that, up to pos- 
sible interference effects, the X(3872) line shapes in the 
near-threshold region can be well described by simple Fla- 
tte formulae, like Eq. ([3]). 
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A The pole structure and the w-plane 
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Fig. 5. Riemann sheets in tj-plane. 



The structure of the poles in the two-channel case can be 
visualised with the mapping of the four-sheeted Riemann 
surface onto the single complex plane |16j . For a given 
energy E one has: 



E 



k\ 



2/V 



E 



k 2 

2^2 



A, 



(58) 



so that, instead of two complex momenta fci and k 2 con- 
strained by two conditions (|58p . it is convenient to intro- 
duce one complex variable lj defined such that 



fci 



HiA 



uj + - 

U! 



1 

LJ 



(59) 



which gives for the energy: 



A ( a 1 
E = — (oj 2 + — 

4 V LJ 2 



By construction, the complex w-plane is free of unitary 
cuts. 

In Fig. [S]we show the mapping of Riemann sheets onto 
the uj complex plane. Thick solid line corresponds to the 
real values of the energy E on the first sheet. The sheets 
arc labelled as follows: 



I 

II 
III 

IV 



Im fci > 0, Im k 2 > 0, 

Im fci < 0, Im k 2 > 0, 

Im fci > 0, Im k 2 < 0, 

Im fci < 0, Im k 2 < 0. 



